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ABSTRACT: The complete coding sequence of Wild Argali ISG15 cDNA was generated by rapid amplification of cDNA ends. The 
ISG15 cDNA was 642 bp with an open reading frame of 474 bp, which encoded a 17.47 kDa protein composed of 157 amino acids. Its 
amino acid sequence shared 97.9%, 80.8%, 91.4%, 94.3%, 78.3% identity with those of ISGlScDNA from Ovis aries (accession no. 
NMOO 1009735.1), Capra hircus (accession no. HQ329186.1), Bos taurus (accession no. BC102318.1), Bubalus hubalis (accession no. 
HM543269.1), and Sus scrofa (accession no. EU647216.1), respectively. The entire coding sequence was inserted into the pET-28a 
vector and expressed in E. coli. The recombinant protein corresponded to the expected molecular mass of 25 kDa as judged by SDS- 
PAGE, and it was detected in the bacterial inclusion bodies. The expressed protein could be purified by Ni^^ chelate affinity 
chromatography and the results from the lymphocyte proliferation test showed that the product could stimulate lymphocyte proliferation 
very well (p<0.05), which further confirmed its biological activity. (Key Words: Argali, ISG15cDNA, Clone, Sequencing, Expression, 
Activity Detection) 



INTRODUCTION 

Argali {Ovis ammon), also known as the big sheep, 
Artiodactyla, Ruminantia, Boridae, Caprinae, genus Ovis, 
belongs to wild animals under the second class state 
protection (Yu et al., 2008). Wild Argali are distributed in 
the mountains of Eurasia, and also found in Inner Mongolia, 
Ningxia, Gansu, Xinjiang, Qinghai, and Tibet in China. 
Their main characteristics are: the rams weigh up to 200 kg, 
both male and female have horns, the whole body is brown 
gray with a white hip spot, they have more lean and less fat 
content enabling fast growth, strong adaptability and high 
disease resistance as (Jueken et al., 2007). 

One of the most abundantly expressed proteins in 
response to interferons (IFN) is IFN- stimulated gene 
product 15 (ISG15), which is not only seen during antiviral 
challenge but also during early pregnancy within the uterus 
of bovine, ovine, porcine, murine, human and baboon 
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species (Rempel et al., 2005). ISG15 is induced as a disease 
resistance gene and expressed efficiently in mammalian 
cells after stimuli with IFNs and viruses (Der et al., 1998). 
ISG15 is a 17.15 kDa protein composed of 165 amino acid 
residues with further modifications after translation. ISG15 
is related to disease resistance in sheep. As congenital 
immune response is the first line of defense against 
pathogens mainly through a series of specific cytokine - 
induced expression (Chang et al., 2006), the rapid high 
expression of ISG15 protein indicates that it may play an 
important role in innate immunity during infection. 

However, there are no reports regarding the ISO 15 
sequence of Xinjiang Wild Argali. The encoded gene 
product corresponds to ubiquitin (Zhao et al., 2004), and the 
amino acids of ubiquitin are not conserved among different 
species (Zhang et al., 2011). To study the difference in 
nucleotide sequences and amino acid sequences of ISG15 
gene between wild and domestic animals, the present study 
used the sequences from cattle and sheep as the references 
to clone and analyze the cDNA sequence by RACE 
technique and in vitro expression in prokaryotic cells, which 
provides a basis for further studies on the differences in 
gene functions and disease resistance. 
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MATEREALS AND METHODS 

Peripheral blood materials and RNA preparation 

Heparinised blood (2 mg/mL) was collected by 
venipuncture from four healthy wild Argali. The peripheral 
blood lymphocyte cells (PBLCs) were separated by using 
the standard FicoU-hypaque gradient centrifugation method 
(Haoyang Biological Manufacturer Co., LTD, Tianjin, 
China) according to the manufacturer's instructions. Total 
RNA was extracted by using the Trizol (Invitrogen, Beijing, 
China) method and stored at -70°C until use. 

Rapid amplification of cDNAS'ends (5'-RACE) 

The 5 ' ends were generated according to the instructions 
for the 5' Full RACE Core Set cDNA Kit (Clontech, Beijing, 
China). Briefly, cDNA was synthesized by reverse 
transcription (RT) UPM primer (as provided in the kit). 
Based on the published mRNA sequence of Ovis aries 
(accession no. NMOO 1009735.1), Capra hircus (accession 
no. HQ329 186.1) and Bos taurus (accession no. 
BC102318.1), gene specific primers namely Bl 
(5TTTGCTACAACTTTATTACTG30, B2 (5'CTTTTCAA 
AAG-ACAGCCAGA3'), CI (5'AAGCAGTGGTATCAAC 
GCAGAGTACGCGGGG3'), and C2 (5'-CCACTCAAA 
CTCCTGCACCAGCCTACC3') were employed in the 
nested PCR at concentration of 20 pmol/|iL each. The 
second amplification products were recovered using Takara 
Agarose gel DNA purification Kit version 2.0' according to 
the manufacturer's instructions and ligated into pMD-lST 
cloning vector (Takara Biotechnology, Dalian, China) as 
described (Sambrook et al., 2002). Seven clones with 
inserts of the expected size were identified after EcoRI and 
Hindlll digestion. The sequences in the seven clones were 
verified by the dideoxy-mediated DNA sequencing method 
using Taq DNA Polymerase (Qiagen, Shanghai, China). 

Rapid amplification of cDNA 3'ends (3'-RACE) 

The 3'ends were synthesized according to the 
instructions in 3 'Full RACE Core Set cDNA Kit (Takara 
Biotechnology). The 3 -sites adapter Primer (provided in the 
Kit), the D3 primer (5'GGGACC-TGAAGGTGAAG 
ATGCT30 and D4 primer (5TGCCTCTGAGGGA 
CTCCATGAT3') were used at a concentration of 20 
pmol/|uL each. The 3 ' end PCR products were purified and 
ligated into pMD-18T cloning vector (Takara 
Biotechnology). Six clones with inserts of the expected size 
were identified after EcoRI and Hindlll digestion. The 
sequences in the six clones were verified by the dideoxy- 
mediated DNA sequencing method using Taq DNA 
Polymerase (Qiagen, Shanghai, China). 

Synthesis of complete ISG15 cDNA 

The complete coding sequence was generated by RT- 



PCR as described (Sambrook et al., 2002). The primers 
used were as followed: PI (5'GAATTCATGGGCG 
GGGACCTGA3'); P2 (5'AAGCTTTTTG-CTACAACTTT 
ATTCACTGCG3') at concentrations of 20 pmol/|iL each. 
The purified PCR products were ligated into pMD-18T 
cloning vector (Takara Biotechnology, Dalian, China), 
transformed and sub-cultured. Ten randomly selected clones 
with inserts of the expected size were identified after EcoRI 
and Hindlll digestion. The sequences of the ten clones were 
also verified by the dideoxy-mediated DNA sequencing 
method using Taq DNA Polymerase (Qiagen, Shanghai, 
China). 

Construction of expression vector 

The expression vector of recombinant ISG15 protein 
was constructed by inserting full length ISG15 cDNA into 
pET-28a plasmid (Takara Biotechnology, Dalian, China) as 
described (Sambrook et al., 2002). Briefly, recombinant 
pMD-18T with ISG15 cDNA was digested with EcoRI and 
Hindlll and the ISG15 cDNA fragment was inserted into 
pET-28a vector previously linearized with similar enzymes 
to obtain the pET-28a/ISG15 cDNA expression vector. The 
recombinant vector was transformed into Escherichia coli 
and positive clones were verified by enzyme digestion. 

Expression of recombinant ISG15 protein 

Transformed Escherichia coli harboring pET-28a/ISG15 
cDNA was sub-cultured in Lurria Bertani media 
supplemented with Kanamycin (100 |ig/mL) and incubated 
at 37°C until the ODgoo reached 0.4 to 0.6. The expression 
was induced with isopropyl-B-D-thiogalactopyranoside 
(IPTG) (Sigma, Shanghai, China) at a final concentration of 
1 mM. After 8 h of incubation at 37°C, the cells were 
harvested and the expression level was assayed by SDS- 
PAGE. 

Purification of recombinant ISG15 proteins 

Following induction, bacterial pellets were collected 
and recombinant proteins were isolated under denatured 
conditions as described by the manufacturer (Qiagen, 
Shanghai, China). Briefly, the bacterial pellet was mixed in 
8 M urea, 0.1 M Na2HP04, 0.01 M Tris-HCl (pH 8.0), 1 M 
2-mercaptoethanol, sonicated and centrifuged. The 
supernatant was applied to Ni-NTA agarose. The NTA 
colunm was washed with 20 volumes of 8 M urea, 0.1 M 
Na2HP04 and 10 mM Tris-HCl (pH 6.3), followed by a 20 
column volume of wash step with 8 M urea and 50 mM Tris 
(pH 8.0). Recombinant ISG15 proteins were refolded in 0.1 
M urea, 50 mM Tris-HCl, 2 M oxidized glutathione (Sigma, 
Shanghai, China), 0.02 M reduced glutathione (Sigma, 
Shanghai, China) pH 8.0. Protein renaturation was carried 
out on the Ni-NTA colunm using a linear gradient over a 
period of 2 h. After renaturation, the proteins were eluted 
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with 0.3 M imidazole in the refolding buffer. The eluted 
proteins were first dialyzed to remove the glutathione in 0.1 
M urea and 50 mM Tris-HCl (pH 8.0), followed by dialysis 
in 0.05 mM urea and 50 mM Tris-HCl (pH 8.0). Protein 
concentration was determined as previously described 
(Bradford, 1976). 

Determination of the activity of the expression product 
by MTT 

Lymphocytes were isolated from the blood of Xinjiang 
Wild Argali using the separation medium (TBD, Ltd, 
Tianjing, China) after incubation for 24 h. The cells were 
inoculated on the 96- well plates at 1x10^ cells/well, 20, 40 
and 80 |ug/mL purified products were added, respectively. 
Meanwhile, another two 96-well plates were set as the 
control group, (Con A was added as the negative control 
(Solarbio Biotechnology, Ltd, Shanghai), while nothing was 
added the blank control). The three 96-well plates were in 
37°C, 5% CO2 for 24 h, then 5 mg/mL MTT (Amresco 
company, Shanghai, China) was added in each well for 
another 4 h. Supernatant was removed and the precipitated 
formazan was dissolved with 100 |uL dimethyl sulfoxide 
(Wolsen company, Xi'an, China) followed by gentle 
shaking of the dishes. Optical density of the dissolved 
formazan was measured at 570 nm by using Thermo 
Labsy stems (Roche, Shanghai, China). Finally, statistical 
analysis was carried out for the data by using one-way 
ANOVAby SPPSS13.0. 

RESULTS 

RACE and complete ISG15 cDNA PGR products 

Nucleotide sequencing combined with a search of the 
GenBank database revealed five clones generated by 
5'RACE-PCR (^500) (Figure 1), as having no the initiator 
element in their open reading frames. The 3'RACE-PCR 
products (Figure 2) had complete 3 'ends containing a 
polyadenylation tail of 12A's at 67 bp downstream of the 
stop codon TAG. All the ten transformation clones had 

1 2 



\< — oOObp 



Figure 1. The 5-RACE-PCR products. Lane 1, DNA maker I 
molecular marker; Lane 2, positive 3 'ends PGR products «500 bp. 
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Figure 2. The 3-RACE-PCR products. Lane 1, DNA maker I 
molecular marker; Lane 2, positive 5 'ends PGR product »560 bp. 

inserts of the expected ISG15 cDNA-ORF as identified by 
EcoRI and Hindlll enzyme digestion (Figure 3). The 
sequences were confirmed in all the ten clones. The 3', 
5 'RACE and full length sequences were 560 bp (3 'ends), 
500 bp (5'ends) and 642 bp (ISG15 cDNA) respectively 
(data not shown), each having 99% identity to Ovis aries 
mRNA. 

Characteristics of ISG15 cDNA sequence 

The sequence of Wild Argali ISG15 cDNA and the 
inferred amino acids (Figure 4) were submitted to GenBank 
as the accession no. KC481400. The accuracy of the 
sequence was determined on both strands and found to be in 
order. The cDNA insert was 642 bp long. ISG15 cDNA- 
ORF contains one long open reading frame of 474 bp from 
90 nt to 563 nt. The first methionine codon ATG begins at 
90-92 nt and ends at presumed stop codon TAG at 561-563 

1 2 3 

Figure 3. The recombinant plasmids of five clones containing 
inserts of the full-length ISG15cDNA-0RF identified by EcoRI 
and Hindlll enzyme digestion. Lane 1, DNA maker I molecular 
marker; Lane 2-3, «541 bp ISG15-0RF. 
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1 AC ATGGGGGAGC GC T GAGAGC GC TGAGAGC GGAGAGGC GGC T AC AGCC AAC C AGT GC C TG 
6 1 CAGAGAACC TGGCACCAGAACCCACGGCCATGGGCGGGGACC TGAAGGTGAAGATGC TAG 
20 MGGDLKVKML 
121 GGGGC GAGGAGAT C C T GGT GC C T C T GAGGGAC TCCATGATGGTATCC GAGC T GAAGC AGT 
40 G GE EI LV PLRDSMM VSE L KQ 
181 TCATCGCCCAGAAGATCAATGTGCCTGCTTTCCAGCAGCGCCTGGCCCACCTTGACAGCA 
60 F lA QK INVP AFQQ RL AH LDS 

2 41 GGGAGGT GC T GC AGGAAGGGGT GC C C C T T GT C C AC C AGGGC C T GAAAGC T GGC AGC AC C A 
80 R EV LQE G VPL VHQ GL K AGST 
301 TCCTGCT GAT GGT GC AGAAC T GC AGC GC C AC C C T GAAC AT CC T GGT GAGGAAC GAC AAGG 
100 ILLMV QN CS AT LNI L VRND K 
361 GT C GC AGCAGC T C C T ATGAGGT C C AGC TGAC GC AGAC T GTGGC T GT GC T C AAGC AGC AGG 
120 GRSS SYE VQ LT QTV AVL KQ Q 

4 21 T GT GC C AAAAGGAGC GTGT GC AAGC GGAC C AGT T C T GGC TGTCTTTT GAAGGGAAGC CCA 
140 V C QK E RVQ ADQ FW LS FEGKP 
481 T GGAT GATGAGC AC C C GC T GGGGGAAT AC GGC C T C AC GACC GGGT GC AC C GT GT T C AT GA 
160 MD DEHP LG FY GLT TG CTVFM 

5 41 ATCTACGTCTGCGGGGTGGGTAGGGAAGGGCCAGGAGGGCCTTAGGGAGGGCTCCCCATG 
180 NLRLR GG * 

601 CAGCGTAGTGAATAAAGTTGTAGCAAAGCCAAAAAAAAAAAA 

Figure 4. Nucleotide and predicted amino acid sequence of ISGlScDNA. The complete ISGlScDNA sequence consisting of 642 
nucleotides and open reading frame of 474 bp. The initiation ATG and stop TAG codons are shaded. The distal polyadenylation signal 
sequence AATAAA are shaded. The nucleotide and its amino acids are numbered along the left margins respectively. Conserved residues 
of the sequence at G, V, M, I, DL, QK, VKML, VSELK, AFQQRLA are shaded. 



nt. The bases AATAAA at 611-626 nt downstream are the 
distal polyadenylation stop signals for ISG15 cDNA. 
However, the TATA box was not found in the upstream of 
ISG15 cDNA -ORE 

The initiation ATG and stop TAG codons are shaded. 
The distal polyadenylation signal sequence AATAAA are 
shaded. The nucleotide and its amino acids are numbered 
along the left margins respectively. Conserved residues of 
the sequence at G, V, M, I, DL, QK, VKML, VSELK, 
AFQQRLA are shaded. 
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Figure 5. Expression of ISG15 protein. Gels were stained with 
Coomassie blue and 10 |liL samples were loaded perlane. Lane 1, 
pET-28a empty expression vector (negative control); Lane 2, E. 
coli (BL21 strain) extracts transformed into pET-28a, before 
induction with IPTG; Lane 3: standard protein molecular weight 
marker; Lanes 4, 5, 6, 7, 8, 9, 10, 11, and 12, transformed E. coli 
(BL 21strain) extract by pET-28a, after 0.5, 1, 1.5, 2, 2.5, 3, 4, 5, 6 
h of induction with IPTG, respectively. 



Predicted amino acid sequence of recombinant ISG15 
cDNA 

The open reading frame encoded a 157 amino acid 
residue protein with an isoelectric point (pi) of 7.242 and a 
predicted molecular mass of 17,471.33 Da (^17.47 kDa). 
The recombinant protein corresponded to the expected 
molecular mass of 25 kDa as judged by SDS-PAGE and the 
expression product was detected in the bacterial inclusion 
bodies (Figures 5 and 6). The protein consists of 54 
hydrophobic amino acids (Ala, lie. Leu, Phe, Trp, and Val), 
40 Polar amino acids (N, C, Q, S, T, Y) (Asn, Cys, Gin, Ser, 
Thr, and Tyr), 17 strongly basic amino acids (Lys and Arg) 
and 17 strongly acidic amino acids (Asp and Glu). 
Alignment of amino acid sequences of ISG15 cDNA against 

1 2 3 



31kDa- 
20.1kDa- 



-25kDa 



Figure 6. Purification of ISG15 protein. Gels were stained with 
Coomassie blue and 10 |liL samples were loaded per lane. Lane 1, 
standard protein molecular weight marker; Lanes 2 and 3, purified 
recombinant 25 kDa, ISG15 protein. 
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Figure 7. Computer genetic distance from ISGlScDNA sequence. Sequences were derived from Gen bank databases. 



amino acid sequences from other species revealed ISO 15 
protein shared 97.9%, 80.8%, 91.4%, 94.3%, 78.3% 
identity with those of ISG15 cDNA from Ovis aries 
(accession no. NM001009735.1), Capra hircus (accession 
no. HQ329186.1), Bos taurus (accession no. BC102318.1), 
Bubalus bubalis (accession no. HM543269.1), and Sus 
scrofa (accession no. EU647216.1), respectively (Figure 7). 
The phylogenetic tree generated from nucleotide sequences 
(Figure 8) revealed a close relationship between Argali and 
Ovis aries. 



Table 1. Effects of Argali 


ISG15 protein on lymphocyte 


proliferation (X+S), (p<0.05) 




Concentration (|Lig/mL) 


OD57O nm 
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20 (ISG15) 


0.1410+0.0057"^ 


40 (ISG15) 


0.1763+0.0061^ 
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Effects of recombination protein on lymphocyte 
proliferation 

The results from lymphocyte proliferation test showed 
that the expression product could stimulate lymphocyte 
proliferation very well (p<0.05), which proved the 
expressed protein had biological activity (Table 1). 

DISCUSSION 

The ISG15 cDNA sequence of the Xinjiang Wild Argali 
for testing was obtained by RT-PCR, 5 'RACE, and 3 'RACE 
methods. The ISG15 cDNA was 642 bp long with an open 
reading frame of 474 bp, which encoded a 17.47 kDa 
protein composed of 157 amino acids, containing one 
initiation codon ATG and one termination codon TAG. The 



5'UTR ends at 89 nt and the 3'UTR ends at 79 bp in ArgaH, 
while the 5'UTR ends at 91 nt, the 3'UTR ends at 81 bp in 
sheep. We found that Wild Argali, goats and sheep showed 
the highest homology by BLAST, which was in accordance 
with our expectation. The bases AATAAA at 611-626 nt 
downstream are the distal polyadenylation stop signals for 
ISO 15 cDNA. However, the TATA box was not found in the 
upstream of ISG15 cDNA-ORF. Song reported the results 
for ISO 15 cDNA cloning from small tail cold sheep and 
jining grey goats where the 3' UTR retained CdxA binding 
sites for transcription even in the case of frame shift, CdxA 
components were generally located in the 5' end of the gene 
and worked as a part of the promoter (Song, 2007). All of 
these results suggest that the CdxA transcription factor may 
regulate the expression of ISG15 gene. The difference in the 
non-coding regions of ISO 15 gene between wild argali and 
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Figure 8. Computer generated phylogenetic tree from ISG15cDNA sequences. Sequences were derived from Gen bank databases. 
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sheep suggests that there may be differences in the 
regulation of ISG15 expression by CdxA transcription 
factor among different species. 

ISG15, as a disease resistance related gene, plays an 
important role in regulating innate immunity (Okumura et 
al., 2006). Ritchie et al. confirmed that ISG15 and covalent 
modifications in ISG15 played important roles in 
stimulating lymphocyte proliferation, enhancing cellular 
antiviral action and promoting leukocyte chemotaxis 
(Ritchie et al., 2004). Moreover, ISG15 is able to induce 
IFN-y of T cell expression, stimulate proliferation of natural 
killer cells, and improve its cytotoxic effect (Kim et al., 
2003). ISO 15 can regulate immune cell activation and 
proliferation (Pitha-Rowe et al., 2007). Therefore, MTT 
method was used to detect the activity of the recombinant 
ISG15 protein in this study. 

The results showed that the recombinant ISG15 protein 
can stimulate lymphocyte proliferation and had certain 
biological activity. Previously, in order to make full use of 
Xinjiang wild argali's superior performance and the 
characteristics of high disease resistance, argali were 
crossed with Xinjiang local sheep. The hybrids grew fast 
and had a larger and better production performance. 
However, it was found that hybrids suffered from poor 
disease resistance, they were susceptible to respiratory 
disease, and the mortality was high during early growth and 
development. 

Cui found after comparing the ISG15 gene sequences of 
parents' generation and hybrids that there were differences 
among them (Cui, 2012). Therefore, the present study was 
carried out to clone and express the ISO 15 gene of the 
argali, which will provide a basis for further studies on the 
differences in gene functions and disease resistance. 
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